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Abstract. In the first part of this series we characterized all linear opera- 
tors on spaces of multivariate polynomials preserving the property of being 
non- vanishing in products of open circular domains. For such sets this com- 
pletes the multivariate generalization of the classification program initiated by 
Polya-Schur for univariate real polynomials. We build on these classification 
theorems to develop here a theory of multivariate stable polynomials. Appli- 
cations and examples show that this theory provides a natural framework for 
dealing in a uniform way with Lee- Yang type problems in statistical mechan- 
ics, combinatorics, and geometric function theory in one or several variables. 
In particular, we answer a question of Hinkkanen on multivariate apolarity. 
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Introduction 

In two seminal papers from 1952 [30|,l61j Lee and Yang proposed the program of 
analyzing phase transitions in terms of zeros of the partition function and proved a 
celebrated theorem locating the zeros of the partition function of the ferromagnetic 
Ising model on the imaginary axis in the complex magnetic plane. This theorem has 
since been proved and generalized in many ways by e.g. Asano, Fisher, Newman, 
Ruelle, Lieb-Sokal, Biskup et al, etc; see ^jSland references therein. Nevertheless, the 
Lee- Yang theorem seems to have retained an aura of mystique. In his 1988 Gibbs 
lecture [48J Ruelle proclaimed: "I have called this beautiful result a failure because, 
while it has important applications in physics, it remains at this time isolated in 
mathematics." Ruelle's statement was apparently motivated by the fact that the 
Lee- Yang theorem also inspired speculations about possible statistical mechanics 
models underlying the zeros of Riemann or Selberg zeta functions and the Weil 
conjectures [551 [351111] but "the miracle has not happened" 

Recently, Lee- Yang like problems and techniques have appeared in various math- 
ematical contexts such as combinatorics, complex analysis, matrix theory and prob- 
ability theory [11 [6l H E [9. JO, 13, M,, 16, 21, M, 25, 41, ■ The past decade has 
also been marked by important developments on other aspects of phase transitions, 
conformal invariance, percolation theory |271 1291 [55] . However, as Hinkkanen noted 
in [25] , the power in the ideas behind the Lee- Yang theorem has not yet been fully 
exploited: "It seems that the theory of polynomials, linear in each variable, that 
do not have zeros in a given multidisk or a more general set, has a long way to go, 
and has so far unnoticed connections to various other concepts in mathematics." 

In this paper we show that the Lee- Yang theorem and the mathematics around it 
are intimately connected with the dynamics of zero loci of multivariate polynomials 
under linear transformations and Problems [T]^ below. As we point out in SjHl such 
connections have been implicitly noted in essentially all known proofs and exten- 
sions of the Lee- Yang theorem. For instance, Lieb and Sokal f37 reduced the at the 
time best Lee- Yang theorem, due to Newman [3H], to the following statement: if 
P,Q G C[zi, . . . , Zn] are polynomials which are non- vanishing when all variables are 
in the open right half-plane, then the polynomial P{d/dzi, . . . , d/dz„)Q{zi, . . . , Zn) 
also has this property unless it is identically zero (Theorem 18. 3|) . Thus, to better 
understand Lee- Yang type theorems one is naturally led to consider the problems 
of describing linear operators on polynomial spaces that preserve the property of 
being non- vanishing when the variables are in prescribed subsets of C". 

Let us formulate these problems explicitly as in 5J. Given an integer n > 1 and 

C C" we say that / S C[2:i, . . . , z„] is fl-stable if /(zi, . . . , Zn) ^ whenever 
(zi, . . . , Zn) G O. A K-linear operator T : V ^ IK[zi, . . . , z„], where K = R or 
C and is a subspace of . . . ,Zn], is said to preserve VL- stability if for any 

ri-stable polynomial f & V the polynomial T(/) is either 17-stable or T{f) = 0. For 
K = e N" let Kk[zi, . . . ,z„] = {/ G K[zi,...,z„] : deg^^(/) < < 

i < n}, where deg^. (/) is the degree of / in z^. By slight abuse of terminology, if 
^ C C and = '5" then il-stable polynomials will also be referred to as ^f-stable. 

Problem 1. Characterize all linear operators T : M^k[zi, ■ ■ ■ , Zn] — > . . . , z„] 

that preserve ^l-stability for a given set C C" and n G N". 

Problem 2. Characterize all linear operators T : K.[zi, . . . , Zn] K.[zi, . . . , Zn] 
that preserve fl-stability, where Q. are prescribed subsets o/C". 
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In physics [541 156j it is useful to distinguish between hard-core pair interactions 
(subject to constraints, e.g. the maximum degree of a graph) and soft-core pair 
interactions (essentially constraint-free) . By analogy with this dichotomy, one may 
say that results pertaining to Problem[l]are "hard" or "algebraic" (bounded degree) 
while those for Problem[2]are "soft" or "transcendental" (unbounded degree), cf. [5]. 

For n = 1, K = R, and = {z e C : lm{z) > 0} Problems [THS amount 
to classifying linear operators that preserve the set of real polynomials with all 
real zeros. This question has a distinguished history that goes back to Hermite, 
Laguerre, Hurwitz and Polya-Schur, see ^ and references therein. In particular, in 
|41j Polya and Schur characterized all diagonal operators with this property, which 
led to a rich subsequent literature on this subject [6l[8llT7 t lT8 l [T9 l[20l[2Tl[26l l3T 1 
l35l l40l l43l [52] , However, it was not until very recently that full solutions to this 
question - and, more generally, to Problems [iHl] for n — I and any open circular 
domain 17 - were obtained in 6J . Quite recently. Problems [THU were solved in [S] 
whenever fl = fli x ■ ■ ■ x il„ and the fi^'s are open circular domains. For such sets 
these results complete the multivariate generalization of the classification program 
initiated by Polya-Schur [41]. They also go beyond e.g. [6l[8] and have interesting 
consequences, as we will now see. 

In Part A we build on the classification theorems of to develop a self-contained 
theory of multivariate stable polynomials. To begin with, in fJT] we study opera- 
tors on multi-afhnc polynomials inspired by natural time evolutions (symmetric 
exclusion processes) for interacting particle systems [34]. We give a new simple 
proof of [ini Theorem 4.20] (see also [33]) stating that these operators preserve 
stability and extend it to all circular domains. In !j2] we use these symmetriza- 
tion procedures to give a new proof of the Grace- Walsh-Szego coincidence theorem 
that unlike most proofs known so far avoids (univariate) apolarity theory. In [J3] 
we establish a "master composition theorem" that extends to several variables all 
the classical Hadamard-Schur convolution results due to Schur-Malo-Szego, Walsh, 
Cohn-Egervary-Szego, de Bruijn, etc [ITl [35l [43] . This also generalizes the multi- 
variate composition theorems based on the Weyl product f5^ as well as Hinkkanen's 
theorem 25 and provides a unifying framework for results of this type. In Sj4]we ob- 
tain "hard" multivariate generalizations of Polya-Schur's classification of multiplier 
sequences that extend the "soft" theorems of (8l. 

As noted in ^44j , the concept of apolarity has a rich pedigree going all the way 
back to ApoUonius and was much studied in invariant theory, umbral calculus, and 
algebraic geometry [55!, ^531. In [l¥j Rota adds: "Grace's [apolarity] theorem is an 
instance of what might be called a sturdy theorem. For almost one hundred years 
it has resisted all attempts at generalization. Almost all known results about the 
distribution of zeros of [univariate] polynomials in the complex plane are corollaries 
of Grace's theorem." In fj5] we establish Grace type theorems for multivariate 
polynomials and provide an answer to a question of Hinkkanen [251 §5]. In ^ 
we prove "hard" Lieb-Sokal lemmas that sharpen the "soft" ones in [32] (whose 
importance in the Lee- Yang program is explained in §8.11 ) 

In Part B we study statistical mechanical and combinatorial applications of the 
theory of stable polynomials developed in Part A and [5]. We show that the key 
steps in existing proofs and extensions of the Lee- Yang and Heilmann-Lieb theorems 
as well as various other theorems on graph polynomials follow in a simple and vmified 
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manner from this theory. These results are due to Asano [2] , Ruelle [45l |46l [5T] , 
Newman [37l|3g, Lieb-Sokal [32], Hinkkanen 25J, Choe et al [16|, Wagner [59] . 



Let us first fix some of the notation that we will use throughout. Recall that the 
support of a polynomial f{z) = X^qgN" G C[zi, . . . , z„] is the set supp(/) = 

{a e N" : a{a) ^ 0}, where z = {zi, . . . , z„) e C", a = (ai, . . . , a„) G N", and 

= zj*! • ■ • z,^". Set [n] = {1, . . . , n} and (1") = (1, . . . , 1) G N". We say that / is 
of degree at most k, = (ki, . . . , k„) G N" if / G Ck[zi, . . . , z„] (cf. the introduction) 
and of degree k if deg^.(/) = k^, i G [n]. Polynomials in C(iii)[zi, . . . , z„] are called 
multi-ajfine. 

We employ the usual partial order on N": if a = (ai, . . . , a„) G N" and (3 = 
. . . , /3„) G N" then a < /3 if < /3i for all i G [n]. Let \a\ = ai + . . . + a„, 
a! = ai! • • • a„!, and 



The open unit disk is denoted by D and open half-planes bordering on the origin 
by He = {z G C : Im(e*^z) > 0}, where 9 G [0, 27r). Note that Mq is the open upper 
half-plane while Mjl is the open right half-plane. Ho-stable polynomials are referred 
to as stable polynomials and those with all real coefficients are called real stable^ 
cf. [5]-[T^. We denote the sets of stable, respectively real stable polynomials in 
n variables by W„(C), respectively Hn(IR)- Polynomials in C[zi, . . . , z„] which are 
H.I -stable are said to be weakly Hurwitz stable. In |16| these are termed polynomials 
with the half-plane property. The notions of Hg-stability are equivalent modulo 
rotations for complex polynomials but this is not so for real polynomials. However, 
for real polynomials with non-negative coefficients PU*, Theorem 4.5] yields the 
following hierarchy of half-plane properties: if such a polynomial is Ho-stable then 
it is He-stable for any 6* G [0, tt]. 



The symmetric group on n elements, S„, acts on C[zi, . . . , z„] by permuting the 
variables: o-(/)(zi, . . . , z„) = , . . . , z^(„)), a G 6„, / G C[zi, . . . , z„]. Define 

the symmetrization operator Sym : C[zi, . . . , z„] C[zi, . . . , Zn] by 



Clearly, Sym is a linear operator whose image consists of symmetric polynomials, 
that is, polynomials invariant under the action of S„. 

Recall that a circular domain in C is any open or closed disk, exterior of a disk, 
or half-plane. The Grace- Walsh-Szego coincidence theorem is an important and 
very useful result on the geometry of polynomials, see, e.g., [TB ] [22 ] [35 l l43 l [60] . 

Theorem 1.1 (Grace- Walsh-Szego). Let f G C(iii)[2i, . . . ,Zn] be a symmetric poly- 
nomial, C be an open or closed circular domain, and ^i, . . . G C. Suppose fur- 
ther that either the total degree of f equals n or that C is convex (or both). Then 
there exists at least one point ^ G C such that 



A. Theory of Multivariate C-Stable Polynomials 




1. Symmetrization Procedures 



Sym(/) = - E ^(/)- 



/(a,...,^n) = /(^,...,c)- 



(1.1) 
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This theorem was essential for proving the sufficiency part of the characterization 
of linear operators preserving C-stability in [5]. We will see here that Theorem 1 1.11 
is actually a consequence of stronger (asymmetric) symmetrization procedures on 
stable polynomials which were used in [lOj to prove correlation inequalities for 
symmetric exclusion processes. More precisely, we will deduce Theorem 11.11 from 
the following result. 

Theorem 1.2. Let C be an open or closed circular domain. 

(a) // C is convex then the symmetrization operator Sym preserves C-stability 
on multi-ajfine polynomials, i.e., Sym : C(i.i)[zi, . . . , z„] — + C(in)[zi, . . . , z„] 
preserves C-stability. 

(b) // C is non-convex and f G C(i»i) [zi, . . . , z„] is C -stable and such that all 
variables are active in f (i.e., df /dzi ^ 0, i E [n]) then Sym(/) is C-stable. 

Remark 1.1. If C is non-convex and / G C(iii) [zi, . . . , z„] is C-stable then the 
condition in Theorem 11.21 (b) that all variables are active is actually equivalent to 
the requirement that / has total degree n, i.e., d^f/dzi ■ ■ ■ 9z„ 7^ 0. 

Remark 1.2. One can easily construct examples showing that Sym does not preserve 
C-stability when acting on arbitrary (not multi-affine) polynomials. 

It is not difhcult to prove Theorem 1 1.2 1 assuming the Grace- Walsh-Szego theorem. 
However, in ^we will prove the latter via Theorem 1 1.2 1 This will make the theory 
developed here and in [5] self-contained. 

We will derive Theorem 11.21 from the next result which was first proved in [TOl 
Theorem 4.20]. From a physical viewpoint [31], Proposition 11.31 implies that stabil- 
ity is preserved by the natural time evolution of symmetric exclusion processes. 

Proposition 1.3. Let f G C(iii) [zi, . . . , z„], C be an open or closed circular do- 
main, < p < 1, and T = (ij) G &„ be a transposition. 

(a) If C is convex and f is C-stable then so is pf + (1 — p)T(f). 

(b) If C is non-convex and f is C-stable and depending on both Zi and Zj then 
p/ + (1 — p)T{f) is also C-stable. 

Our proof of Proposition ! 1 . 31 relies on the maximum principle for harmonic func- 
tions which we use to prove the following lemma. Another recent elementary proof 
of Proposition 11.31 was independently given in ^33] . Let Hq be the closed upper 
half-plane. 

Lemma 1.4. Let f{z, w) = a + bz + cw + dzw G C[z, w] and define 
Vi{f){x) = Im(ac) + Im(aJ + bc)x + lm{bd)x'^ , 
V2{f){x) = Im(a6) + Im(aJ + cb)x + Im(cJ)a;^. 

Suppose that c? 7^ 0. 

(1) /// ts Mo-stable then lm{b/d) > or Im(c/rf) > 0. 

(2) /// IS Wo-stable and Im(c/d) > then Vi{f){x) > and V2(f){x) > for 
all x eR. 

(3) /// IS Wo-stable and lm{b/d) > then V^2(/)(a;) > and Viif){x) > for 
all X gR. _ 

(4) Iflni[b/d) > and Im(c/d) > then f is Mo-stable if and only if for some 
(and then any) i G {1,2} one has Vi{f){x) > for all a; G M. 
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Proof. Since the partial derivative of a Ho-stable polynomial is Ho-stable or identi- 
cally zero (if in doubt apply Theorem l3.ip we have lm{c/d) > and lm{b/d) > if / 
is Ho-stable. Iflm(c/d) = lm{b/d) = then the polynomial (/"^/(z—c/d, w — fe/c?) = 
wz + {ad — bc)/d'^ is Ho-stable. This is a contradiction since {wz : w,z £ Ho} = C. 
Hence Im(c/d) > or lm{b/d) > if / is Ho-stable. 

Assume that lm{b/d) > and Im(c/rf) > 0. Solving for w in f{z,w) = we see 
that / is Ho-stable if and only if 



Now, p is a harmonic function in the half-plane {z G C : Im(z) > — Im(c/d)} which 
contains Ho. Let Kr = {z e C : Im(z) > 0, \z\ < r}. By the maximum principle 
the minimum of p on Kr is attained on the boundary of Kr- For real x we have 



Moreover, p{z) lm{b/d) > as z — > oo. Since the same arguments apply if one 
instead solves for z in f{z,w) = 0, this verifies (4). 

If just Im(c/ d) > we will still have p{x) > for all x G R and if only Im(6/ d) > 
then p{x) > for all x G IR\ {— c/d}. By symmetry in z and w this verifies (2) and 



Recall the multivariate version of Hurwitz' theorem on the "continuity of zeros" , 
of. [M Footnote 3, p. 96]. 

Theorem 1.5 (Hurwitz' theorem). Let D be a domain (open connected set) in C" 
and suppose {fk}kLi « sequence of non-vanishing analytic functions on D that 
converge to f uniformly on compact subsets of D. Then f is either non-vanishing 
on D or else identically zero. 

To deal with discs and exteriors of discs we also need Lemmas 11.61 and ll .81 below 
- which were proved in [5', Lemmas 6.1 and 6.2] - and Corollarv ll.7l 

Lemma 1.6. Let {C'i}^^i be a family of circular domains, f G C[zi, . . . , z„] be of 

degree k G N", and J C [n] a (possibly empty) set such that Cj is the exterior of 
a disk whenever j G J. Denote by g be the polynomial in the variables Zj, j G J , 
obtained by setting Zi = Ci £ Ci arbitrarily for i ^ J. If f is Ci x ■ ■ ■ x Cn-stable 
then swpp(g) has a unique maximal element 7 with respect to the standard partial 
order on N'^ . Moreover, 7 is the same for all choices of Ci Ci, i ^ J . 

An immediate consequence of Lemma ll.6l is the following. 

Corollary 1.7. Let k G N" and 1^ '. Ck[zi, . . . , z„] — > Ck[zi, . . . , z„] be the linear 
operator defined by I^^z") — z""", a < k. Then /„ restricts to a bijection between 
the set of ^-stable (D-stable) polynomials in Ck[zi, . . . , z„] and the set ofC\D-stable 
(C \ D-stable) polynomials in Ck[zi, . . . , z„] of degree k. 

If {Ci}"^]^ is a family of circular domains and k — (ki, . . . , k„) G N" we let 
7Vk(Ci, . . . ,Cn) be the set of Ci x • ■ • x C„-stable polynomials in Ck[zi, . . . , z„] that 
have degree kj in zj whenever Cj is non-convex. Note that if all Cj are convex then 
A/'k(Ci, . . . , Cn) consists of all Ci x • • • x C„-stable polynomials in Ck[zi, . . . , z„]. 




(1.2) 



p{x) 



Im(ac) + lm{ad -|- bc)x + lm{bd)x^ 



\c-\-dx\'^ 



(3) of the lemma. 



□ 
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Recall that a Mobius transformation is a bijective conformal map of the extended 
complex plane C given by 

0(C) = -^q—^, a,b,c,d £ C,ad — be — I. (1-3) 

Note that one usually has the weaker requirement ad — be =/= but this is equivalent 
to (|1.3p which proves to be more convenient. 

Lemma 1.8. Suppose that Ci, . . . , C„, Di, . . . , D„ are open eireular domains and 
K = . . . , Kn) G N". Then there are Mobius transformations 

CjC + di 

as in (|1.3p sitc/i t/iai i/ie invertible ) linear transformation : [zi , . . . , z„] ^ 
Ck[zi, ■ • • , Zn] defined by 

$^(/)(zi, . . . ,Z„) = (ciZi • • • (C„2„ +d„)''"/(0i(zi), . . .,(j)n{Zn)) 

restricts to a bisection between 7Vk(Ci, . . . , C„) and Nk{Di, . . . , Z3„). 



Proof of Proposition \1.3[ Clearly, it is enough to prove the proposition for closed 
circular domains. Suppose first that C is the closed upper half-plane Hq and let 
/ G C(i-n)[zi, . . . ,Zn] be Ho-stable. Assuming, as we may, that i — 1 and j = 2, we 
need to prove that 

Pfi^l , 6, . . . , 60 + (1 - 6, ■ • ■ , in) + 

whenever 6, . . . , ^„ G Hq. By fixing ^3, . . . , arbitrarily in Hq and considering the 
multi-affine polynomial in variables z\ , Z2 given by 

{zx.zi) ^g{zi,Z2) -.^ f{zi,Z2, £,3, 

we see that the problem reduces to proving that for any p € (0, 1) the polynomial 
pf{z,w) + (1 —p)f{w,z) is Ho-stable provided that f{z,w) is Ho-stable. This is 
easy to check if d = so we may assume that d ^ 0. Now, if {i,j} = {1, 2} then 



V,, 



{pf{z, w) + {l- p)f{w, z)) = py, (/(z, + (1 - P)V, [f{z, w)) 



which proves the proposition for C = Ho by Lemma [1.41 (since p G (0,1) implies 
that we will be in case (4) of Lemma 11.41) . 

Let C be a closed disk and suppose that / is C-stable. Then by compactness / is 
C-stable for some open disk C* D C. The result now follows by applying Lemma [TTSl 
(with K = (1"), Ci = C, = D, £ € [n], where D is an arbitrarily fixed open 
half-plane) and using the fact the partial symmetrization operator commutes with 
the operator defined in Lemma [L8l 

The case of the closed exterior of a disk follows from the disk case considered 
above and Corollary [TT7] for k = (1") (cf. Remark ll.ip . □ 

In the theory of interacting particle systems [34j it is well known that the sym- 
metrization of a polynomial / can be achieved by applying />—>(/ -f t(/))/2 
infinitely many times with different transpositions r. For the sake of completeness, 
we will give a proof of this fact in the Appendix. 
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Proof of Theorem \1.2[ Consider first the case when C is an open circular domain 
and recall Remark [TTTl Since Sym(/) is obtained by applying f ^ [f + t(/))/2 
infinitely many times with different transpositions r (see Lemma 19.21 below) the 
result follows from Hurwitz' Theorem 11.51 and Proposition 1 1.31 
If C is closed write 

d 

f{z,...,z) = B\{{z-c^), 
where Cj ^ C for j e [d\ and B ^ 0. Clearly, the polynomial . . . , z„) defined 

by 

d 

F{zi,...,Zn) = Bl[{z,-c,) 

is D-stable, where D is a suitable open circular domain containing C but none of 
the Cj's. Then by the above Sym(i^) is ZJ-stable and since Sym(F) = Sym(/) the 
theorem follows. □ 

2. The Grace- Walsh-Szego Coincidence Theorem 

Using Theorem ll.2l we can give a new proof of the Grace- Walsh-Szego coincidence 
theorem that does not rely on apolarity theory as do most known proofs so far, see 
[22l [35l [43l [57l [60] and SjS) We actually prove a more general version of this result 
which holds for families of circular domains. 

Theorem 2.1. Suppose /(^n, . . . , ^iki, . . . , ^ni, • • ■ , 2nK„) is a multi-affine polyno- 
mial in \k\ complex variables which is symmetric in {zij : j G ['ti]} for all i G [n], 
where n — . . . , k„) G N" with > 1, i G [n]. Let further Ci, i G [n], be circular 
domains and ^.y G C^, j G [Ki], i E [n]. Then there exist S,i £ Ci, i £ [n], such that 

/(Cll I • ■ ■ I I ■ ■ • ? Cnl J • ■ ■ I ^iiK„ ) — /(Cl : ■ ■ ■ 1 S,l : ■ ■ ■ : £,n ■ ■ ■ : ^n) 

provided that f has total degree Ki in {zij : j G [Ki]} whenever Ci is non-convex. 

Proof. Clearly, it is enough to prove that if the polynomial in n variables 

g(zi, . . . , z„) := f{zi, . . . , zi, . . . , z„ . . . , z„) 

is Ci X • • • X C„-stable then /(zn, . . . , . . . , z„i, . . . , z^^J is Cf' x • • • x Q"- 
stable, which we will now do by considering one variable at a time. By assumption 
g has degree Ki in the variable Zi whenever Ci is non-convex (symmetry prevents 
cancellation). Fix Q G Cj, j G [rt — 1]. The polynomial h{zn) '■= g{Cij ■ • ■ i Cn-ij ^n) 
is C„-stable and we may therefore write 

d 

h{Zn) = B Y\_(Zn - Ctj), 

where B and aj ^ Ci for j G [d] , so the polynomial 

d 

H(wi, . . . , Wk„) := B ]^(wj - aj) 

is also C„-stable. Now, if C„ is non-convex then by Lemma 1 1.61 one has d — k„ and 
by Theorem 11.21 the symmetrization operator Sym acting on k„ variables maps Hg 
to a C„-stable polynomial. Since the numbers (^i, i G [n — 1], were arbitrary this 
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means that the polarization of g that spUts the variable Zn symmetrically into Kn 
new variables, i.e., the linear operator 

g{zi, ...,Zn)^ f{zi, . . . . . . ,2„_1, . . . ,Z„_i,Z„i, . . ..ZukJ 

preserves the stability in question. By polarizing one variable at a time we conclude 
that /(zii, . . . , ziKi, . . . , Zni, . . . , 2;„K„) is ' x • • • x C^"-stable. □ 

Remark 2.1. In [47j Ruelle produced a proof of the Grace- Walsh-Szego coincidence 
theorem (Theorem II. ip using similar ideas. 

Remark 2.2. A yet more general version of Theorem l2.1l was actually given by Walsh 
in [60, Theorem 1] without assuming any degree conditions, the only requirement 
in ^60j being that C^, i G [n], should be just (closed) circular domains. However, in 
such generality Walsh's aforementioned result fails already for rt = 1. 

3. Master Composition Theorems 

Composition (or convolution) theorems such as the Schur-Malo-Szego theorems 
f [T7l Theorem 2.4], [43l Theorem 3.4. Id]), the Cohn-Egervary-Szego theorem ( [43l 
Theorem 3. 4. Id]), Walsh's theorems (]43l Theorems 3.4.2c and 5.3.1]) or de Bruijn's 
theorems |1HI12) play an important role in the analytic theory of univariate complex 
polynomials and allow to locate their zeros in certain circular domains [35l |43] . 

Using results of [5] we establish "master composition theorems" that provide a 
unifying framework for multivariate generalizations of the classical theorems men- 
tioned above. Let us first recall two of the classification theorems from [5]. 

Theorem 3.1. Let k, G N", T : Ck,[zi, . . . , z„] — ^ C[zi, . . . , z„] be a linear operator, 
and C — He for some < 9 < 2t:. Then T preserves C -stability if and only if 

(a) T has range of dimension at most one and is of the form T{f) — a{f)P , 
where a is a linear functional on Ck[zi, . . . , z„] and P is a C -stable poly- 
nomial, or 

(b) The polynomial (in 2n variables) 

r[(z + z.)«]:=^ Qr(z")^«-" (3.1) 

a<K 

is C'-stable. 

Remark 3.1. Theorem 13.11 triviallv implies the following (well-known) multivariate 
Gauss-Lucas theorem: if / G [zi , . . . , z„] is H^-stable for some < 6* < 27r 
then df/dzi is Hg-stable or identically zero for any i G [n]. More generally, the 
{n l)-variable polynomial / + Zn+idf/dzi is Hg-stablc. 

Theorem 3.2. Let K.eW\T :CJ zi, . . . , Zn] — > C[zi, . . . , Zn] be a linear operator, 
and C = D or . Then T preserves C -stability if and only if 

(a) T has range of dimension at most one and is of the form T{f) = a{f)P , 
where a is a linear functional on Ck[zi, . . ■ , Zn] and P is a C -stable poly- 
nomial, or 

(b) The polynomial (in 2n variables) 

r[(l + zH1:-E (3-2) 

is C -stable. 
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The polynomials in (13. ip and (|3.2p are called the algebraic symbols of T with 
respect to the circular domains under consideration (for Hz. it is often more con- 
venient - but equivalent - to choose (|3.2|) rather than (|3.ip . cf. 5, Remark 6.1]). 

The main result of this section is as follows. 

Theorem 3.3. Let k = (ki, . . . , Kn) £ N", f{u, v) e C[iii, . . . , u„, vi, . . . , Vn], and 
g{z,w) G C[zi, . . . , z„, wi, . . . ,Wn]- Suppose that deg„. (/) < Hi and deg^.{g) < Ki 
for all i £ [n]. 

(a) // / and g are Mg-stable for some < 6 < lix , then the polynomial (in An 
variables) 

is Wg -stable or identically zero. 

(b) // / and g are Mq- stable, then the polynomial 

T^i-^r'-^-gJ^M-^iz,..) 

is Mfj-stable or identically zero. 

(c) // / and g are O-stable, then the polynomial 

a.<K 

is O-stable or identically zero. 

Proof. We only prove (a) since the proofs of (b) and (c) are almost identical. Let 
7 = (71, ... , 7„) G N" be fixed. Define a C[ui, . . . , w„]-valued linear operator T on 
the space of all polynomials h in 4ri variables wi, . . . , w„ satisfying deg„^ {h) < kj, 
deg„^ {h) < 7j , deg^^. (h) < 7^ , and deg^^. {h) < 7^ for all j G [n] by setting 

T[h{u,v,z,w)] = • d^^^^^'"^^- 

Let u = (t2i, . . . ,Un), V = (ui, . . . ,Vn), z = (zi, . . . ,z„) and w) = (wi, . . . ,Wn) be 
new sets of variables. The symbol of T with respect to Hg, that is, 



= + v)'^{z + z)^ {w + wy g{z + u + u,w) 
is clearly Hg-stable which proves (a) by Theorem l3 . ll since 7 G N" was arbitrary. □ 

An important special case of the above theorem is particularly attractive, as is 
its proof. 

Corollary 3.4. Let k G N" and /, g G C[zi, . . . , 2:„, wi, . . . , w„] be of the form 
fiz,w) = Y 9{z,w) = E 
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(a) // / and g are We -stable for some < < 27r, then so is the polynomial 

unless it is identically zero. 

(b) // / and g are Mo-stable, then so is the polynomial 



unless it is identically zero. 
(c) // / and g are H-stable, then so is the polynomial 

E(:)..wo.w4e^-£((o.^)-0(^>o) 

unless it is identically zero. 

Proof. Suppose that /, g are as in part (a) of the corollary. Let 

[zi , . . . , z„] and S : [zi , . . . , z„] — > C-^ [zi , . . . , z„] 

be the linear operators whose algebraic symbols with respect to Hg (cf. (|3.ip ) are /, 
respectively g, with /3,7 G N" appropriately chosen. By Theorem 13 . 1 1 both 5 and T 
preserve Hg-stability, hence so does their (operator) composition ST whose symbol 
is precisely the polynomial in (a). Applying Theorem 13.11 again we conclude that 
this polynomial is Hg-stable unless it is of the form A{z)B{w) for some polynomials 
A and B. If this is the case and these polynomials are not identically zero then 
A{z) must be Hg-stable (being the polynomial P in part (a) of Theorem 13. 1|1 and 
by exchanging the roles of / and g we get that B{w), thus also A{z)B{w), must be 
Hg-stable. This proves (a). Parts (b) and (c) follow similarly. □ 

Example 1. Let us show how the classical (univariate) Schur-Malo-Szego theorem 
can be easily derived from CoroUarv 13.41 If J22=o iTl^kz'' and X]fc=o iTj^'^z'' are 
two polynomials with real zeros only and in addition all the zeros of the latter 
polynomial are non-positive then the bivariate polynomials 

fc=0 ^ k=0 ^ ^ 

are stable (note that / actually depends only on z). Corollarv l3.4l (a) implies that 
the Schur-Malo-Szego composition of the two given polynomials, i.e., the univariate 
polynomial 

V ( ? lafc&fez'' 



k=0 ^ 

is also stable (that is, real-rooted). 

Example 2. Theorems 3.11 and 4.6 in 8^ (the former actually follows from the 
latter) provide some multivariate extensions of the classical composition results 
mentioned above. In particular, |8j Theorem 4.6] shows that the Weyl product 
of polynomials (defined via the product formula in the Weyl algebra) preserves 
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stability, that is, if f{z,w) and g{z,w) are stable polynomials in z (zi, . . . ,z„) 
and w = (wi, . . . , Wn) then so is the polynomial 

aeN" 

To see that this is in fact a consequence of Theorem 13.31 (b) let kjv = (iV, . . . ,N) £ 
N" and let / and g be stable polynomials as in the statement of Theorem 13.31 
Then since stability is closed under scaling the variables with positive numbers the 
polynomial 

H^iu, V, z,w) = Y: (-1) • g^iNu, V) ■ —iz, w) 

is stable for large N . But then the polynomial 

(_\\a Q<^ j d"g 
Yy^^H^iu/N, v,z,w)^Y. 5^1^^"' ' d^^^'"^^ 

a. 

is stable or identically zero, as claimed. 

We conclude with one further consequence. For t ~ (ii, . . . , i„) G M" define the 
t- deformed Weyl product of / and g by 

Using again the fact that stability is closed under scaling the variables with posi- 
tive constants we deduce from above that the t-deformed Weyl product preserves 
stability. As a special (univariate) case, note that if n = 1, t < and / S ]R[z] \ {0}, 
g e M[i(;] \ {0} have all real zeros then by the above the univariate polynomial 

has all real zeros. We thus recover de Bruijn's 11, Theorem 2] and [l2l Lemma 1]. 

4. Hard Polya-Schur Theory: Bounded Degree Multiplier Sequences 

Using the characterization of linear operators preserving real stability obtained in 
[5] we can establish "hard" (bounded degree) multivariate versions of Polya-Schur's 
classification of multiplier sequences |41j that extend the "soft" (unbounded degree) 
theorems of [B] . 

A sequence of real numbers {A(A;)}fegN is called a multiplier sequence if the linear 
operator on univariate polynomials defined by T{z'') — \[k)z^, fc G N, preserves 
real-rootedness, that is, T{f) G 7ii(K) U {0} whenever / G 7ii(K). A multivariate 
multiplier sequence is then defined as a sequence {A(a)}Q,gN" of real numbers such 
that the linear operator T : R[zi, . . . , 2;„] R[zi,...,z„] defined by T{z°') = 
A(a)z", a G N", preserves real stability, see [5]. These were characterized in [5] but 
here we will prove the corresponding "hard" theorems. Given k G N" we say that a 
sequence {X{ol)} a<K of real numbers is a K-multiplier sequence if the linear operator 
T : Rk[zi, . . . , Zn] ^ Rk[zi, . . . , Zn] defined by T{z") = A(a)z", a < k, preserves 
real stability. This is the multivariate generalization of n-multiplier sequences |19) . 

Recall the following lemma from [S]. 
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Lemma 4.1. Let f,gE M[zi, . . . , z„] \ {0}, set h ~ f + ig and suppose that f and 
g are not constant multiples of each other. The following are equivalent: 

(1) h is stable; 

(2) \h{z)\ > \h{z)\ for all z — (zi, . . . , z„) G C" with Ini(zj) > 0, j G [n], where 

z — , . . . , Z^ ) J 

(3) / + z„+i5GW„+i(C); 

(4) / and g are stable and 

whenever z = (zi, . . . , z„) G C" with lm{zj) > 0, G [n]. 

A polynomial / G C[zi, . . . , z„] \ {0} is said to have the same-phase property 
if there exists a G K such that all the non-zero coefficients of e~*"/ are positive. 
The next lemma was first proved in [16i Theorem 6.f]. We will use it in the proof 
of Lemma 14.31 below and for completeness we provide here a self-contained proof 
based on our results so far. 

Lemma 4.2. // / G C[zi^ . . . , z„] \ {0} is homogeneous of degree d and Mg-stable 
for some 9 G [0, 27r) then f has the same-phase property. 

Proof. Note first that the assumptions of the lemma imply that / is He/ -stable for 
any 6' G [0, 2tt). Without loss of generality we may also assume that djf ^ 0, where 
dj — ^ . We will now use induction on d. The statement is trivially true for d ~ 

so suppose d > 1. Applying ^ to the identity f{tzi, . . . , tz„) = t'^f{zi, . . . , Zn) and 
setting t = 1 we get 

n 

f{zi, . . . ,z„) = d^^^Zjdjf{zi, . . . ,z„). (4.1) 

Each polynomial djf is stable (e.g. by Theorem 13. ip and homogeneous of degree 
d — 1. By the induction hypothesis there exists aj G M such that e^^^Wjf has 
all non-negative coefficients. In view of (|4.ip it is therefore enough to show that 
aj = ttfe mod 27r, j, fc G [n]. For each j G [n] we get by Remark 1 3. 11 Lemma l4Tl (3) 
<^ (4) and homogeneity that 

Im(^Mf)^<0^ z = (zi,...,z„)GH^, (4.2) 

Re(^^^)>0, z = (zi,...,z„)eH|. (4.3) 

By continuity (|4.2p also holds for all z G M" for which f{z) ^ 0. Using homogeneity 
we see that lm{dj f {— z) / f {— zj) = —lm{djf{z)/f{z)) hence djf{z)/f{z) is a real 
rational function. Since e~'^"'djf{(^) G M we deduce that e~'"j /(C) £ R for j G [n] 
and C e K"- Now, from gj]) with z G R!^ and the fact that e'^'^Wjfiz) > for 
all such z we conclude that e~'"j/(z) > whenever z G M", j G [n], and thus 
ttj = Qffe mod 27r, j, fc G [n], as required. □ 

Lemma 4.3. Let f{z, w) — X^asN" o,{a)z"^w°' G C[zi, . . . , z„, Wi, . . . , w„] . T/ien / 
is stable if and only if it can be written as 

f{z, w) = C'fi{ziwi) ■ ■ ■ /„(z„w„), 



and 
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where C G C and fi{t), . . . , fn{t) are univariate real polynomials with real and 
non-negative zeros only. 

Proof. The sufficiency part follows simply by noticing that if /i < then /i + zw is 
a stable polynomial in two variables. 

Suppose that /(z, w) = X^qsN" o,{a)z°'w°' is stable. We claim that its support 
J := supp(/) has unique minimal and maximal elements with respect to the stan- 
dard partial order on N". Assume the contrary, let a, a' be two different minimal 
elements and let i, j be indices such that > a'^ and aj < a^. By 13, Theorem 3.2] 
J is a jump system and since a" = a — ei ^ J, there is an index k ^ i such that 
a" + ek J. Let g{zi, Wi, Zk, Wk) be the polynomial 

dw"" 

and set 



ze=wi=0, £^{i,k} 



h{zi,Wi,Zk,Wk) = lim A g{\zi, Xwi, \zk, Xwk). 

By Hurwitz' theorem (Theorem II .51) ft, is a stable polynomial. However, by con- 
struction h is of the form AziWi + BzkWk with AB ^ 0, which is a contradiction 
since polynomials of this type cannot be stable. This shows that J has a unique 
minimal element. 

If f{z, 1) has degree at most Ki in the variable Zi, i e [n], we may consider the sta- 
ble polynomial z'^w'^ f{—z^^ , —w^^), where k — (ki, . . . , k„), z~^ — , . . . , z~^) 
and similarly for w^^. By the above the support of the latter polynomial has a 
unique minimal element, thus providing a unique maximal element for the support 
Jof /. 

Let now ^, k be the minimal, respectively maximal element of J and let T : 
Ck[zi, . . . , Zn] —>■ Ck[zi, • . . , Zn] be the linear operator defined by 

T[z"] = A(a)z" = (-1)" Q a(a)z", < a < k. 

Let {ej}"^i be the standard orthonormal basis of R". We want to show that 

A(a) =A(e)-"+iA(C+(ai-a)ei)---A(e + K-en)e„), ^ < a < k. (4.4) 

This will then prove the necessity part since / will split into a product as in the 
statement of the lemma and the polynomials fj{t), I < j < n, will have the desired 
properties since fj^ZjWj) is necessarily stable. 

To prove (|4.4p note that the algebraic symbol of T is given by 

Gt{z,w) = ^(-l)"a(a)z"w''-" = w'^/(z,-u;-i), 

which is stable. By Theorem 13.11 T preserves stability. Since Gt{z,w) is homoge- 
neous we may assume that (— l)"a(a) > for all a in view of Lemma 14.21 Now, it 
is easy to check that 

a + bz + cw + dzw € M.[z, w] is stable <;=^ be > ad (4-5) 

(see, e.g., [13] or just adapt the arguments in the proof of Lemma fO)) and of course 

a + 2bz + cz^ G K[z] is stable 6^ > ac. (4.6) 

Let 7 G N" and 1 < i,j < n be such that j + Ci + ej < k. Applying T to the 
polynomials 2^(1 -I- Zi){l + Zj) and z^{l — Zi){l + zj) and keeping (|4.5p and (|4.6p in 
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mind we see that A(7)A(7 + ei + Cj) > \{j + ei)A(7 + Cj) and A(7)A(7 + + Cj) < 
A(7 + ei)A(7 + Cj) hence 

Ml)Hl + ei + ej) = A(7 + ei)A(7 + ej) whenever 7 e N" and -y + e. + Cj < n. (4.7) 

From gJl) and 113^ CoroUary 3.7] we deduce that A(7) > for all ^ < 7 < k. The 
proposed formula ()4.4p now follows by induction over k |a| — □ 

Recall the following theorem from 5 . 

Theorem 4.4. Let k e N" and T : Mk[zi, . . . , z„] — > M[zi,...,z„] be a linear 
operator. Then T preserves real stability if and only if either 

(a) T has at most 2 -dimensional range and is given by T{f) = a{f)P + f3{f)Q, 
where a, (3 are real linear forms on Rk[-Z1j • ■ • , Zn] and P,Q€ H„(R) are 
such that P + iQ £ Hni'C), or 

(b) Either T[{z + wY] G 7^2«(K) or T[{z - w)"] G 7^2«(K). 

The "hard" multivariate version of Polya-Schur's theorem [IT] is as follows. 

Theorem 4.5. Let k G N", let A :— {X{a)}a<K be a sequence of real numbers 
and T : Mk[zi, . . . , z„] Mk[zi, . . . ,Zn\ be the corresponding linear operator. The 
following are equivalent: 

(a) X is a ^-multiplier sequence; 

(b) ±A is the product of one- dimensional Ki-multiplier sequences that are either 
all alternating in sign or all non-negative, i.e., 

±A(q;) = Ai(ai) • • • A„(a„), a < k, 

where \i is a Ki-multiplier sequence, i G [n], and either all Xi's are non- 
negative or all are alternating in sign; 

(c) Etther T[{z + w^] G 7^2«(K) or T[{z - wY] G H2„(M); 

(d) T[{z + wY] or T[{z — wY] can be written as 

fl{ziWi) ■ ■ ■ fn{ZnWn), 

where fi{t), . . . , fn{t) are univariate polynomials with real zeros only, and 
all these zeros have the same sign (collectively). 

Proof. This is an immediate consequence of Theorem 14.41 and Lemma 14.31 □ 

Remark 4.1. Note that the k- multiplier sequences with constant sign are precisely 
the sequences whose corresponding operators preserve stability. 

Corollary 4.6. Let k G N", A :— {X{a)}a<K be a sequence of complex numbers 
and T : [zi, . . . , Zn] Ck [zi, . . . , Zn] be the corresponding linear operator. Then 
T preserves weak Hurwitz stability if and only if X is (a constant complex multiple 
of)a non-negative K-multiplier sequence. 

Proof. By Theorem 13.11 T preserves weak Hurwitz stability if and only if the poly- 
nomial 

T[{z + wY] = J2 (^)Ha)z"w''-" 
is weakly Hurwitz stable, which occurs exactly when 

a<K ^ ^ 

is stable. The assertion now follows from Lemma l473l and Remark 1 4. II □ 
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5. Multivariate Apolarity 

The goal of this section is to develop a higher-dimensional apolarity theory and 
establish Grace type theorems for arbitrary multivariate polynomials. 

Two univariate polynomials f{z) = J2k=o (fc)'^fc-^'^ ^^"^ •9(^) ~ Y^k=Q ik)^f^^'' 
degree at most n are apolar if 

{f:9}n E(-i)V<'=^(o)g("-^-)(o) = ;^E(-i)'Ur''^"-^- ^ ^- 

*;=0 ■ fe=0 ^ ^ 

Grace's classical apolarity theorem is as follows 1551 H51 [57] . 

Theorem 5.1 (Grace). Let f and g be apolar polynomials of degree n > 1. If f 
has all zeros in a circular domain C then g has at least one zero in C. 

Note that we may reformulate Grace's theorem as follows. 

Theorem 5.2. Let f and g be polynomials of degree n > I and C be a circular 
domain. If f is C-stable and g is <C \ C-stable then {f,g}n ^ 0. 

For two polynomials /, g £ C[zi, . . . , z„] and k G N" define 

{f,9}. :-E(-ir/(")(0).9('=-"H0) 

and call / and g apolar if they both have degree at most k and {/, g}^ = 0. 

Hinkkanen [25] wondered if Grace's theorem could be extended to several vari- 
ables (he actually only considered multi-affine polynomials) but the precise form 
of such an extension remained uncertain. He also claimed that arguments due to 
Ruelle and Dyson [171 110] could be extended to prove the following result. 

Lemma 5.3. Let A and B be closed subsets of <C which do not contain the origin 
and let f,g G C(i2)[zi, Z2]. // / is A x B -stable and g is {<C\A) x {<C \ B)-stable 
then {f,g}(i2) ^ 0. 

Lemma [5.31 is false, as one can see by considering for instance /(zi, Z2) = zi + Z2, 
g{zi, Z2) — 1, and A = B = {Im(2;) > 1}. However, it holds under additional degree 
constraints (e.g. if both / ang g have total degree 2) which are tacitly assumed in 
[50t Footnote 7]. In [23 Hinkkanen also proposed two possible generalizations of 
Grace's theorem as the following questions. 

Question 1 (Hinkkanen). Let Ai, i € [n], be closed subsets of C that do not 
contain the origin and f,g£ C(in) [zi, . . . , z„]. If f is Ai x ■ ■ ■ x An-stable and g is 
(C \ Ai) X • • • X (C \ An)-stable then {/, ^ 0. 

Question 2 (Hinkkanen). Let Ci, i G [n], be closed circular domains and f,g £ 
C(in)[zi, . . . , Zn\. If f is Ci X ■ ■ ■ X Cn-stable and g is (C\Ci) x • • • x {<C\Cn)- stable 
then ^ 0. 

We will now see that these questions are not true in full generality, but if we 
strengthen the hypothesis slightly in the second question then it is true for arbitrary 
degree polynomials. 

Note that if /, g G Ck[zi, . . . , z„] then 

{f,g}.{z) :^^(-l)"/(")(z)g(«-")(z) 

OL<K 
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is a constant function so in this case {f,g}K{z) = {/,<?}« for z 6 C". Elementary 
computations also yield the following. 

Lemma 5.4. Let k = . . . , k„) G N", a^, 6^, G C, a^dj — 6,jCj = 1, i G [n], 
/, g G Ck[zi, . . . , z„] and se^ 

F(z) = (cizi +di) ' • • • (c„z„ + dn) "/ — -, . . . , — — , 

\CiZi+di CnZn + dnJ 

G[z) = (cizi + di) 1 • • • (c„z„ + d„) — — , . . . , — — . 

\CiZi+di CnZn+dnJ 

Then{f,g}^^{F,G},. 

Remark 5.1. Note that Lemma (5.41 asserts that the functional {•,•}«; is invariant 
under the action of the group of Mobius transformations normalized as in (|1.3p . For 
n = 1 this is quite well-known [43j and motivates the name "apolar invariant" for 
{•, - j^ which is classically used in invariant theory, umbral calculus, and the theory 
of algebraic curves [2H1 Sll [53] . 

Lemma 5.5. Let f,g G C[zi, . . . , z„] and suppose that g has degree k G N". // / 
is H-stable and g is C\ D-sia6/e then {/, g}^ ^ 0. 

Proof. Let 

f{z) = ^ a^z" and g{z) = ^ 6„z". 

a a 

Then h{z) := D-stable by Corollarv ll.7l By compactness there 

is e > such that \h{z)\ > e for z G D". This means that there is (5 > such that 
Ea(-l)"^K-a(l + (5)l"lz" is D-stable. Then by applying Corollary [33] (c) to the 
D-stable polynomials 

^a„z"w" and + 

a a 

we deduce that the polynomial 

is D-stable or identically zero. The assumptions on / and g guarantee that aobn ^ 0, 
so F is D-stable and 

{/, g}, = KlF{{l + 5)-''\ . . . , (1-1- <5)-i/2) ^ 0, 

as claimed. □ 

We find it most natural to state two apolarity theorems: one for discs and exterior 
of discs (Theorem 15. 6p and one for half-planes fTheorem l5.8p . 

Theorem 5.6. Let Ci, i G [n], be open discs or exterior of discs and let f,g & 
Ck[zi, . . . , Zn]- Suppose that 

(i) f is Ci X ■ ■ ■ X Cn-stable and deg^^. (/) — Kj whenever Cj is the exterior of 
a disk, and 

(ii) g is {<C\ Ci) x • ■ • x (C \ Cn)-stable and deg^.{f) — Kj whenever Cj is a 
disk. 

Then 7^ 0. 
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Proof. For i G [n] let <i>i be a Mobius transformation {atz + bi) / {az + di) as in ()1.3|) 
for which <i>i(D) = d. By Lemma [1.81 the polynomial 

F(z) := (ciZi + di)"! • • • (C„Z„ + d„)''"/(*l(^l), ■ • ■ , ^n{Zn)) 

is D-stable. Note that —di/ci ^ C \ D since none of the CiS is a half-plane. The 
polynomial 

G{z) := (cizi • •■ (c„z„ + d„)'""5($i(zi),.. .,$„(z„)) 

is therefore C\D-stable. By Lemma [L8l the degree of g is precisely k so Lemma ISTSl 
applies and by Lemma EH] we get {/, = {F, G}k 7^ 0. □ 

The homogeneous part of a polynomial / G C[zi, . . . , z„] is the polynomial fu 
obtained by extracting the terms of maximum total degree, i.e., 

fnizi, ...,Zn)= lim r'^f{tzi, tz^), 

t — ^oo 

where d = max{|a| : (9"//9z" 7^ 0}. 

Lemma 5.7. Suppose J G C[zi, . . . , z„] is Ho-siaWe anrf i G [n]. T/ien 

\dzi) jj dzi ■ 

Proof. Suppose that / has total degree d. Clearly, it is enough to prove that either 
df/dzi is identically zero or its total degree is d — 1. 

Assume that df /dzi ^ and its total degree is rf' < d — 1. By Remark |3. II the 
polynomial / + Zn+idf / dzi is Ho-stable. Consider now the univariate polynomials 

p(z) = /(j,...,z) = /a(l,...,l)z'' + --, 

These polynomials are of degree d and d' respectively, since the leading coefficients 
are non-zero by Lemma 14.21 Solving for z^+i in p{z) -t- Zn+iq{z) = we see that 
1uy{p{z) / q{z)) > whenever Im(z) > 0. This is a contradiction since p{z)/q{z) = 
Cz^-d' ^ o(2<i-rf') when z ^ 00, where C ^ and d - d' > 2. □ 

Theorem 5.8. Let Ci and C2 be two open half-planes with non-empty intersection, 
K G N", and f,g & C„[zi, . . . ,z„]. /// is C'l-stable, g is C2-stable, and k < a -\- P 
for some a G supp(/), /? G supp(g), then {/, g}^ 7^ 0. 

Proof. By an affine transformation of the variables we may assume that there is 
an e > such that f{z — ie) and g{—z + ie) are Ho-stable, where e ~ (e, . . . , e). 
Then so are the 2n-variable polynomials f{z -\- w — ie) and g{—z — w -\- ie) and by 
Corollary 13.41 also the polynomial 

unless it is identically zero. If it is not identically zero then the conclusion of the 
theorem follows by setting z = w = ie. To complete the proof we show that F{z, w) 
is not identically zero. Let Ga(z, w) = - i£)g('=-")(-2 + ie). This 

polynomial is Ho-stable or identically zero (by Remark 13. ip and by Lemma 15.71 
and Lemma 14.21 all non-zero coefficients in its homogeneous part have the same 
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phase as those in the homogeneous part of f{w ~ ie)g{—z + ie). By the as- 
sumptions on the supports of / and g there is an a such that Ga{z,w) 7^ so 
Imit^oot^'''^'^^^^^ F^tz.tw) ^ 0, where d and e are the total degrees of / and g, 
respectively. In particular, F{z, w) is not identically zero. □ 



6. Hard Lieb-Sokal Lemmas 

In [32I Proposition 2.2] Lieb and Sokal proved that the operation that replaces 
one variable with differentiation with respect to another variable preserves weak 
Hurwitz stability. This result played a key role in the study of Laplace transforms 
of Lee- Yang measures and the extensions of Newman's strong Lee- Yang theorem ob- 
tained in "32] , see i J8.1l It was also an essential ingredient in proving the sufficiency 
part of the classification theorems of [5 . 

The Lieb-Sokal result is a "soft" (transcendental/unbounded degree) result since 
it amounts to saying that the linear operator on C-[zi, . . . , Zn] acting on monomials 
as 

^"^(-ir- ^ " a=K,...,a„)eN", (6.1) 

(which one may schematically represent as "zi ^ —d/dz2" ) preserves (]HIo-)stability, 
see Theorem [O] in EH 

By considering certain linear operators on finite-dimensional polynomial spaces 
we can establish "hard" versions of Lieb-Sokal's result. 

Lemma 6.1. Let n,d G N with n > 2 and let k G N" be such that ki = K2 = d. 
Define a linear operator Td : C„[zi, . . . , z„] — > Ck[zi, . . . , z„] by 



Qdf 



dlf-^^dz'^dzf' 
Then Td preserves Wg- stability for any < 9 < 27r. 

Proof. The symbol of Td, i.e., the 2n-variable polynomial Td[{z + w)"] is given by 

i(z3 + Wsr^ ■■■iZn+ WnT- ^ 773fcV + ^if'^'i^^ + ^2)'' 

k=0 ^ ' 

= (Zl -I- Z2 + U;i + W2f{z3 + Ws)"" • ■ ■ (Zn + Wn^" 

which is Hg-stable. The conclusion follows from Theorem 13.11 □ 

Remark 6.1. An interesting property of Td is that Td{f) is actually a polynomial 
in the n — 1 variables zi + Z2, Z3, . . . , z„ for any / G Ck[zi, . . . , z„]. Indeed, let 
F{zi, Z2, . . . , Zn) — Td{f){zi, Z2, . . . , z„). It is straightforward to show that 

d 

— F(zi + t, Z2 - t, Z3, . . . , 2„) 0, i G C, 
at 

and by letting i = Z2 we get F[zi + t, Z2 — t, . . . , Zn) — F{zi + Z2, 0, . . . , 2„). 

Using Lemma 16.11 and Remark 16.11 we deduce the following "hard" result that 
substantially improves (|6.ip when the top degree is specified. 



20 



J. BORCEA AND P. BRANDEN 



Corollary 6.2. Let rt, d e N with n > 2 and let k G N" be such that ki = K2 = d. 
Define linear operators Sd,Rd '■ Ck[zi, . . . , Zn] — > Ck[zi, . . . , z„] by 



and 



' d 
.k=0 



1 , . / a 



k=0 



Rd 



Lfc=0 



= ^h-^nd-ky.{^)\k 

k=0 ^ ^ 



(Z2, . . . ,Z„). 



T/ien S'c; and Rd preserve (Rq-) stability up to degree k. 



The above "hard" results do indeed imply the "soft" ones. To see this fix /3 G N" 
and set (/3)q = for a € N". In [5] Lemma 8.2] it was shown that the linear 

operator on C[zi, . . . , Zn] defined by z" i— > {(3)aZ°', a G W\ preserves stability. 
In particular, if '^f.^Q ZiQk{z2) is stable then X^fc^o (d-fc)! Qfc('^2) is stable and 
this extends to n variables. Therefore, the "soft" Lieb-Sokal result (respectively. 
Theorem I8.3P follows from Lemma 16.11 (respectively. Corollary [ 



7. Transcendental Symbols and the Weyl Algebra 

Define the complex Laguerre-Polya class 7i„(C) as the class of entire functions 
in n variables that are limits, uniformly on compact sets, of polynomials in Ti.n(C), 
see, e.g., [211 Chap. IX]. The usual (real) Laguerre-Polya class 7i„(M) consists of 
all functions in 7i„(C) with real coefficients. 

If T : K[zi, . . . , Zn] — * ]K[zi, . . . , z„], where K = M or C, is a linear operator 
we define its transcendental symbol, Gt{z,w), to be the formal power series in 
wi, . . . , Wn with polynomial coefficients in K[zi, . . . , z„] given by 

Gt{z,w) ^ (-l)"r(z")^. 

oGN" 



By abuse of notation we write Gt{z, w) = T[e~^ '"], where z-w = + . . ■ + ZnWn- 
Let us recall from [5] the transcendental characterizations of complex, respectively 
real stability preservers. 

Theorem 7.1. Let T : C[zi, . . . , z„] — > C[zi, . . . , z„] be a linear operator. Then T 
preserves (Mq-) stability if and only if either 

(a) T has range of dimension at most one and is given by T{f) — a{f)P, where 
a is a linear form on C[zi, . . . , z„] and P G 7i„(C), or 

(b) GT{z,w)&H2n{C). 



Remark 7.1. From Theorem 17.11 one can easily deduce a characterization of linear 
operators preserving fJ-stability for any open half-plane SI. For instance, the analog 
of Theorem 17.11 (b) for the open right half-plane is that the transcendental 
symbol of T with respect to ]HI.|, i.e., the formal power series 

r[e--] := T{z-)^, 

defines an entire function which is the limit, uniformly on compact sets, of weakly 
Hurwitz stable polynomials. 
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Theorem 7.2. Let T : ]R[zi, . . . , 2„] . . . , z„] be a linear operator. Then T 

preserves real stability if and only if either 

(a) T has at most 2 -dimensional range and is given by T{f) — a{f)P + (3(f)Q, 
where a, (3 are real linear forms on M[zi, . . . , z„] and P,Q ^ 7i„(]R) are such 
thatP + iQ G HniC), or 

(b) Either Griz, w) e H2rM) or Gt{z, ~w) e H2„(R). 

To illustrate the power of Theorems 17.11 and 17.21 we show that the main results 
of [5j for partial differential operators immediately follow from these two theorems. 
Recall that a (Weyl algebra) finite order linear partial differential operator with 
polynomial coefficients is an operator T : K.[zi, . . . , 2;„] — * IK[zi, . . . , 2;„], where 
K = C or M, of the form 

a<l3 

where /3 e and Qa G K[zi, . . . , Zn], a < p. 

Theorem 7.3. Let T : &[zi, . . . , z„] ^ '^[zi, . . . ,Zn], where K = C or R, be defined 
by (|7.ip and set 

F[z,w) = ^ Qa{z)w°' G K[2;i, . . . z„, wi, . . . , w„]. 

Then 

(a) T preserves stability if and only if F{z, —w) is stable; 

(b) T preserves real stability if and only if F(z, —w) is real stable. 

Proof. The (transcendental) symbol of T is given by 

T[e-'-'"]^e-'-'"Fiz,-w), 
so (a) and (b) follow immediately from Theorems 17.11 and I7.2[ respectively. □ 

Remark 7.2. Theorem [73] was first established in [51 Theorems 1.2-1.3] by different 
methods. An interesting consequence noted in |8l Theorem 1.11] is that if a Weyl 
algebra operator T preserves (real) stability then so does its Fischer-Fock dual T* . 
As shown in |5] , this duality result is a powerful multivariate generalization of the 
classical Hermite-Poulain- Jensen theorem and Polya's curve theorem [T71 [33] . 

B. Applications 

We will now apply the theory developed in Part A to show that (the key steps 
in) existing proofs and generalizations of the Lee- Yang and Heilmann-Lieb theo- 
rems follow in a simple and unified way from the characterizations of rj-stability 
preservers in terms of operator symbols obtained in [5]. These results are due to 
Asano [5], Ruelle [351 [31 [SI], Newman [371 [35], Lieb-Sokal Hinkkanen [15], 
Choe et al (TH], Wagner [SS]- For brevity's sake, we will only focus on the main 
arguments used in deriving them and in some cases we point out possible extensions. 
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8. Recovering Lee- Yang and Heilmann-Lieb Type Theorems 

Let us first recall the original version of the Lee- Yang theorem for the partition 
function of the ferromagnetic Ising model (at inverse temperature 1). This function 
may be written as 

where a ■ h = X]r=i '^i^i ^^^'^ l^i'^) ~ e^»'j=i -^^i'^^'^i , 

Theorem 8.1 (Lee- Yang [30]). If Jij > for all i.j G [n] then 

(a) Z{hi, . . . , hn) 7^ whenever Re(/ii) > 0, 1 < i < n; 

(b) All zeros of Z{h, . . . ,h) lie on the imaginary axis. 

Remark 8.1. In physical terms [31 SI [301 ISH ISS] j the Jij are ferromagnetic (> 0) 
coupling constants while the hi are external (magnetic) fields sometimes also called 
fugacities. Theorem 18.11 (b) asserts that the zeros of the partition function of the 
ferromagnetic Ising model accumulate on the imaginary axis in the complex fugacity 
plane and a (first-order) phase transition occurs only at zero magnetic field. 

Before we give a proof of the Lee- Yang theorem let us make a historical digression. 
In his work on the zeros of the Riemann zeta function Polya was led to a simple 
yet useful result: 

Lemma 8.2 (Polya [35], Hilfssatz II). Let a > 0, 6 G M, and G{z) he a real entire 
function of genus or 1 with at least one real zero and only real zeros. Then the 
function 

G{z + ia)e'^ + G{z ~ ia)e-''^ (8.1) 

has only real zeros. 

Hilfssatz II was subsequently employed by Kac [40j pp. 424-426] to settle a 
special case of Theorem 18 . 1 1 that proved to be inspirational for Lee and Yang's final 
proof [3D] (cf. [51 Remark 4.2]). Recently, Lee- Yang type results and applications 
to Fourier transforms with all real zeros were obtained in [l] jM] [15] by iterating the 
process of Hilfssatz II. A simple proof of this result and multivariate extensions is 
as follows. Let R be the linear operator on formal power series in n variables with 
complex coefficients f{z) — X^qgN" ol{ol)z" defined by 

i? ( Y «(")^" ) = E Re(a(a))z" ^ i(/(z) +7(1)) . 

\QiEN" / qEN" 

By Lemma 14.11 R maps the set of stable polynomials into the set of real stable 
polynomials and consequently also the complex Laguerre-Polya class Hn (C) (cf . SJ7]) 
into the Laguerre-Polya class 7Y„(M). In the special case when n = 1 and G{z) is as 
in Lemma l8^ it follows from Hadamard's factorization theorem that G{z) £ 7ii(R) 
hence G{z + ia)e^^ £ Hi(C) and by the above 

2R{G{z + ia)e''') = G{z + ia)e''' + G{z - ia)e-''' G Hi(R), 

so the function in (|8.H) has only real zeros. Note also that e^*^ G H„(C) and thus 
2R{e-") = cos(z) G 7?„(R). 

More general versions of Theorem l8 . 1 1 were obtained in e.g. [32] and [38], see §8.11 
For simplicity of argument and exposition we will concentrate for the moment just 
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on the original Lee- Yang theorem and give a short proof based on the ideas in [32] 
combined with Theorem 17. II 

Proof of Theorem \8.1\ Note that (b) foUows from (a) by symmetry in —a. To 
prove (a) define to be the set of functions ji : { — 1,1}" C whose Laplace 
transform 

<Te{-i.i}" 

is the limit, uniformly on compact sets, of weakly Hurwitz stable polynomials (i.e., 
non- vanishing whenever all variables are in the open right half-plane HI.|). 

Claim: Let i,j e [n] and > 0. If /i G then /iy G M.^ where 

e'-''5/i((T) if fjj = o-j. 



e •■''^ /i((t) if (Ti ^ fjj. 



Let us show that the claim implies the theorem. Indeed, if /io : { — 1,1}" — > C 
is such that /x(cr) = 1 for all <j G { — 1, 1}" then its Laplace transform Z^^ equals 
(e'*^ +6^''^) • • • (e''" +e^'*"). As noted above one has cos(z) G 7i„(M), which implies 
that /xq G by a rotation of the variables. Then by successively applying to /io 
the transformations defined above for all pairs («, j) G [n] x [n] one gets (a). 

To prove the claim note that Zj^^^. — T{Z^^), where 

T = cosh(Jij) + smh(J^j) ^^^^ . 

By Theorem 17.11 and Remark 17.11 the operator T preserves weak Hurwitz stability. 
Since T is a second order (linear) differential operator, by standard results in com- 
plex analysis we have that if /fc — *■ / uniformly on compacts then T(fk) — > T{f) 
uniformly on compacts. This proves the claim. □ 

8.1. Newman's Theorem and the Lieb-Sokal Approach. In [38] Newman 
proved a strong Lee- Yang theorem stating that the Lee- Yang property holds for 
one-component ferromagnetic pair interactions if and only if it holds for zero pair 
interactions. This theorem was subsequently generalized in [32] by Lieb and Sokal 
who showed that one-component ferromagnetic pair interactions are "universal mul- 
tipliers for Lee- Yang measures" and established a similar result for two-component 
ferromagnets. Lieb-Sokal's key observation was that it would suffice to show that 
a certain linear differential operator preserves the Lee- Yang property, which they 
proved by reducing the problem to the following statement about polynomials. 

Theorem 8.3 (Lieb-Sokal). Let {Pi{u)}'^i and {Qi{v)}™^i be polynomials in n 
complex variables u = (ui, . . . , Un), v = [vi, . . . , w„), and define 



m 

^(z) = E^^(^/^^)Q'W' 
1=1 

where z ~ (zi, . . . , 2„), d/dz = (d/dzi, . . . , d/dzn). If R is weakly Hurwitz stable 
(in In variables) then S is either weakly Hurwitz stable or identically zero. 
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Proof. Define a linear operator 

T : C[ui, . . . ,Un,Vi, . . . ,Vn] ^ C[ui, ...,Un,Vi,.. .,Vn] 

by letting 

and extending linearly. Clearly, the theorem is equivalent to proving that T pre- 
serves weak Hurwitz stability. By Theorem 17.11 and Remark 17.11 this amounts to 
showing that the formal power series 



a,/3 



a!/3! 



(i.e., the transcendental symbol for H^i) defines an entire function which is the limit, 
uniformly on compact sets, of weakly Hurwitz stable polynomials. An elementary 
computation then yields 



1=1 



which satisfies the above requirement since e^"" — lim„^oc (l + zw/n)^. □ 

8.2. The Schur-Hadamard Product and Convolution. The following version 
of Theorem 18.11 is usually referred to as the Lee- Yang "circle theorem", see, e.g., 
[251 EH. 



Theorem 8.4. Let A = {aij) be a Hermitian n x n matrix whose entries are in 
the closed unit disk D. Then the polynomial 

/(zi, . . . , z„) = n n "-^3 

5C[n] ieSjffS 

is I}-stable. In particular, f{z,...,z) has all its zeros on the unit circle. 

Hinkkanen's proof [25] of Theorem 18.41 makes use of a composition theorem for 
the Schur-Hadamard product of multi-affine polynomials which is defined as follows: 
if /(^) = Escw ci{S)zS and g{z) = Escf„i KS)z'^ then 



(/.g)(z)= <S)h{S)z' 



ayj )U(J 

SC[n] 

The next result is Hinkkanen's composition theorem [25, Theorem C]. 



Theorem 8.5. Let f,g G C(in) [zi, . . . , . If f,g are B)-stable then so is f • g 
unless it is identically zero. 

Proof. Let g be a fixed D-stable multi-affine polynomial in n variables and let T 
be the linear transformation on multi-affine polynomials in n variables given by 
T{f) = f»g. Recall Theorem (b). The symbol of T is 



T 



(l + zu))["l 



.g(ziwi, . . . ,z„w„). 



which is clearly B-stable (in 2n variables). Theorem 13.21 vields the result. □ 
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Remark 8.2. The proof of Theorem 18.41 given in [25j is as follows. For i,j £ [n] 
with i < j (note the typo "i ^ j" in [25]) let 

fij{zi, . . . ,Zn) ^ {1 + aijZi + a[JZj + ZiZj) Y\_ (1 + ^/c)- 

keln]\{i,j} 

It is not hard to see that fij is D-stable and by taking the Schur-Hadamard product 
of all these polynomials one gets 

(/i2 • • • • • /(„-!)„) (2) = XI n n 

SC[n] ieS j^S 

which is again D-stable by Theorem 18.51 



Using Corollary [33] we can extend Hinkkanen's composition theorem to arbitrary 
(not necessarily multi-affine) D-stable polynomials: 

Theorem 8.6. Let f{z) = □a(a)z" and g{z) = OK")^" be B- 

stable polynomials. Then so is 



unless it is identically zero. 

Proof. Apply Corollary [33] (c) to the ©-stable polynomials f{z) and g{zw), where 

ZW = (ziWi, . . . , ZnW„). □ 



Closely related to the Schur-Hadamard product is the convolution operator on 
multi-afhne polynomials jl6j defined as follows: if f{z) = X]sc[n] o,{'S)z^ and giz) 



EscM biS)z' then 



(/*5)(^)- E <S)b{T)z 



SAT 

S,TC[n] 



where S'AT = (5 U T) \ (5 fl T). A corresponding composition result - this time 
for weak Hurwitz stability - is given in [161 Proposition 4.20]. 

Theorem 8.7. Let f,gG C(iti')[zi, . . . , z„]. // f,g are weakly Hurwitz stable then 
so is f -k g unless it is identically zero. 

Proof. Let g be a fixed -stable multi-affine polynomial in n variables and let 
T be the linear transformation on multi-affine polynomials in n variables given by 
Tif) = f*9- The symbol of T (cf. Theorem [33]) is just 



[z + wy^ ^{z + wy^g\ ■ ■ . (8.2) 

J \Zi+Wi Zn + Wn ) 

Now if u, w G Hii then also u''^ , V"^ ^ u + w e hence 

U + V ^ ' ^ 

so that the polynomial in (|8.2p is H:| -stable (in 2n variables). Theorem 13.11 again 
yields the desired conclusion. □ 
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8.3. Asano Contractions. Many known proofs of the Lee- Yang theorem are based 
on so-called Asano contractions or variations thereof [21 HSl [50l [51] . Let 

/(Zl, . . . , Z„) = a(z3, ...,Zn)+ 6(2:3, . . . , Zn)zi + c{z3, . . . , Z„)z2 + d{zz, . . . , Zn)ZiZ2 

be a polynomial in n > 2 variables which is multi-affine in zi and Z2- The Asano 
contraction of / is 

A{f){zi, . . . , z„) = a(z3, . . . ,z„) -I- d(z3, . . . , z„)zi. 

Note that A[f) does not depend on Z2- The key fact used in the aforementioned 
proofs is a property of Asano contractions that may be stated as follows. 

Lemma 8.8. Let k = (ki, . . . , k„) G N" with n > 2 and ki = K2 = 1. Then 

A : Ck[zi, . . . ,z„] Ck[zi, . . . , z„] 

is a linear operator that preserves "O- stability. 

Proof. It is clear that A is linear. Its (algebraic) symbol is 

^[(1 + zwY] = (1 zw)(''^-'''")(l -I- ziu;iW2), 

which is D-stable, so the assertion follows from Theorem 13.21 □ 

8.4. Multi-AfRne Part and Folding mod 2. Recall that a matching in a graph 
G = (y, E) is a subset M of E such that no vertex of the graph {V, M) has degree 
exceeding one. The general version of the Heilmann-Lieb theorem on the monomer- 
dimer model is the following. 

Theorem 8.9 (Heilmann-Lieb [23]). Let G = {V, E) he a loopless graph and define 
its matching polynomial with edge weights {\e\eeE o,nd vertex weights {zijigy as 

AfG(z,A)= ^ ^eZtZj. 

rnatchings M e—ij^M 

If 1^0, e e E, then Mq{z, A) is a weakly Hurwitz stable polynomial (in z). 

In [TB] and [SS] §5] it was shown that the Lee- Yang and Heilmann-Lieb theorems 
can actually be given a unified combinatorial formulation and proof. The idea is to 
form the "test" polynomial 

Fg{z,X)= n i'^ + ^eZrZj) (8.3) 
e={i..j}eE 

associated to a graph G = {V,E), \V\ = n, equipped with vertex weights {zi}i^v 
and non-negative edge weights {AejeS-E- This polynomial is weakly Hurwitz stable 
in the z^'s and by applying to it appropriate linear operators one gets precisely the 
polynomials occurring in the Lee- Yang theorem and the Heilmann-Lieb theorem, 
respectively. Thus one only has to check that the linear operators used in this 
process preserve weak Hurwitz stability. These operators are defined as follows. 

The linear operator MAP : C[zi, . . . , z„] — > C(in) [zi, . . . , z„] extracts the multi- 
affine part of a polynomial, that is, if /(z) = X^aGN" 0(0^)^" then 

MAP(/)(z) = a(a)z". 

a: Qi<l, z£[n] 
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The transcendental symbol (see Remark I7.ip of MAP is 



a: <1, [nl 



Clearly, this is a weakly Hurwitz stable polynomial. Theorem 17.11 and Remark l7. II 
imply that MAP preserves weak Hurwitz stability. It is easy to see that 

Mg(z, A) = MAP \Fg{z, A)j , (8.4) 

which yields the Heilmann-Lieb theorem. 

The linear operator MOD : C[zi, . . . , z„] C(in)[zi, . . . , z„] "folds mod 2" the 
powers in the Taylor expansion of a polynomial, i.e., if f{z) = X^aeN" then 



MOD(/)(z) = «(")^ 



a mod 2 



where a mod 2 = (ai mod 2, . . . , q;„ mod 2). The algebraic symbol of MOD (up to 
degree k) with respect to the open right-half plane (cf. Theorem 13. 2p is 



M0D[(1 + zwT] = 2-"(l + wT{l + z)["l W 



If Re(C) > then |(1 — C)/(l -t- C)| < 1 so the above symbol is weakly Hurwitz 
stable. By Theorem 13.21 we conclude that MOD preserves weak Hurwitz stability. 
MOD is employed in [16l §4.8] to prove the Asano contraction lemma (Lemma l8.8p 
and thereby the Lee- Yang theorem as well (cf. i J8.3l and §5]). 

8.5. Generalizations of the Heilmann-Lieb Theorem. It is natural to study 
graph polynomials with more general degree constraints than those defining the 
matching polynomial and to establish Heilmann-Lieb type theorems for such poly- 
nomials. This has been pursued by e.g. Ruelle [JS] |46j and Wagner [58l|59]. If 
G = {V, E) is a graph with vertex set V = {vi, ...,?;„} we let deg G G be the 
degree vector of G, i.e., the i-th coordinate of degG is the degree of Vi in G. Let 
K — (ki, . . . ,K„) e N" and suppose that degG < k. Then given degree weights 
It : N'' ^ C and non-negative edge weights {)\e]eeE one may ask what are the 
non- vanishing properties of the polynomial 

Fg(z,A,u)= ^ A^u(deg(y,i?))^'^°s(^'^). (8.5) 

H<ZE 

This question was considered by Wagner in [59]. When it comes to weak Hurwitz 
stability it is natural (and of course sufficient) to require that the linear "truncation" 
operator T : Ck[-Zi, . . . , z„] ^ Ci^[zi, . . . , z„] defined by T(z") = u{a)z°' , u{a) G C, 
a < K, preserves weak Hurwitz stability. Now these are precisely the multivariate 
multiplier sequences (up to degree k) that were characterized in Corollary 14.61 of 
this paper as follows: 

u(a) = ui(ai) • • • u„(a„), a < k, 
where for each i S [n] the polynomial 



fe=0 



u,{k)z'' (8.6) 
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has all real non-positive zeros. We thus recover the following generalization of the 
Heilmann-Lieb theorem due to Wagner [59] , which extends a theorem of Ruelle [45] . 
By the necessity in Corollarv l4.6l the theorem below is optimal. 



Theorem 8.10 (Wagner [SH])- Let G — {V,E) be a graph whose degree vector 
satisfies deg G < k and let Fg{z, X,u) and u he as in (jS.Sp and (|8.6p . respectively. 
If {\e}eeE are non-negative edge weights then 

(a) Fq{z, X,u) is weakly Hurwitz stable considered as a polynomial in z; 

(b) All zeros of the univariate polynomial 

E <<^^dv,H))\At^ 

k=0 \ HCE I 

\\H\ = fe / 

are real and non-positive. 

Proof. The first statement follows from Corollary 14.61 and the fact that the test 
polynomial defined in (|8.3p is weakly Hurwitz stable. If we set all the z^ 's, 1 < j < n, 
equal to —it wc obtain the univariate polynomial (in t) 

k=Q \ HQE I 

\\H\ = k / 

which is then real stable. Clearly, this forces the polynomial in (b) to have all zeros 
real and non-positive. □ 

Remark 8.3. In ^59] Wagner also proves non-vanishing properties in sectors, which 
cannot be obtained by our methods. However, Theorem l8.10l is slightly more general 
in that wc consider max-degree at every vertex (not uniform max-degree). 

9. Appendix 

We give here a simple proof of the fact that Sym can be viewed as a (convergent) 
infinite product of operations as those in Proposition 11.31 For a G S„ define an 
operator To- : C[zi, . . . , z„] C[zi, . . . , Zn] by 

TAf) - ^ E 

l\ ''I re{a) 

where (a) is the subgroup of 6„ generated by a. Given a, /3 G N" we write a ^ P if 
a is a rearrangement of /3. If f(z) = X^qsN" o,{a)z°' e . . . , z„] let the symmetry 
index of / be defined by s(/) = X]a~/3 Wi'^) ~ a{(3)\. For f = g -\-ih ^ C[zi, . . . ,Zn] 
with 5, /i G M.[zi, . . . ,Zn] we define its symmetry index as s(/) = s{g)+s{K). Hence 
s(/) = if and only if / is symmetric. 

Lemma 9.1. Let a G 6„ and f{z) = a{a)z°' G C[zi, . . . , z„]. Then 

s(T.(/))<s(/) (9.1) 
with equality if and only if a {a (a)) = a{a) for all a G N", i.e., T„[f) = f. 
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Proof. We may assume that / £ R[zi, . . . ,Zn]- Since s(t(/)) — s{f) for all r 6 6^ 
we have by the triangle inequality 

<TM)) = ^ a(r(a))-a(r(/3)) 

^ E E Hria))-a{rm\ 

E ^(^(/)) 



= <f) 

with equality if and only if the following condition holds: 

If a ~ /9 then a{T{a)) ~ a(r(/3)) have the same sign for all r G (a). (A) 

Clearly, if a{a{a)) — a{a) for all a G N" then equality holds in (|9.ip . On the other 
hand, if equality in (|9.ip holds let /3 = cr(a) and assume that a{a) > a{[3) (the case 
a{a) < a{(3) follows similarly). Then by (A) we have 

a{a) > a{a{a)), a{a{a)) > a{a'^{a)), a{a''~^ (a)) > a{a) , 

where k is the order of a. Hence a a(a) is constant on (CT)-orbits, which completes 
the proof. □ 

For / G C[zi, Zn] let 

A{f) = {Tr^^ ■ ■ -Tr^if) : Ti, . . . ,Tk € &n are transpositions} 

and denote by A{f) the set of polynomials that are limits, uniformly on compact 
sets, of polynomials in A{f). 

Lemma 9.2. If f e C[zi, . . . , z„] then Sym(/) G Aif). 

Proof. We claim that the set s(^(/)) := {s{g) : g G A{f )} is closed. Suppose that 
Xk X as k —>■ oo, where Xk = s(5fc) with G A{f) for fc G N. Let | ■ |r be the 
supremum norm on the ball of radius r in C". If cr G S„ we have by the triangle 
inequality and invariance under permutations that |Tcr((7)|r < \g\r- It follows that 
\h\r < \f\r for all h G A{f). Hence, by Montel's theorem, {gfelfeeN forms a normal 
family so there is a subsequence converging uniformly on compacts to a polynomial 
g G A{f) with 5{g) = x. _ 

Hence y := infs(^(/)) is achieved for some g G A{f). If T{g) ^ g for some 
transposition r G 6n then by Lemma l9. II we have 5(Tr{g)) < s{g). However, one 
clearly has Tr{g) G A{f), which contradicts the minimality of s{g). We deduce that 
T{g) — g for all transpositions r G 6„ and thus g — Sym(/) and 5{g) — 0. □ 
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